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ABSTRACT 


In  a  wedge  shaped  ocean,  the  method  of  images  is  used  to 
develop  an  analytical  approximation  of  the  acoustic  pressure 
field.  Contemporary  work  develops  acoustic  doublets  from  a 
combination  of  the  source  and  surface  reflection  image  using 
simple  dipole  theory.  The  method  of  images  is  then  used  to 
sum  the  dipole  images.  This  thesis  matches  dipole  pairs  to 
achieve  a  quadruplet  expansion.  A  computer  program  using  the 
derived  quadruplet  equation  is  then  created  to  verity  the 
results  by  comparing  them  with  the  "URTEXT"  program. 

DTIC  QUALITY  INSPECTED  8 


iil 


TABLE  OP  CONTENTS 


I.  INTRODUCTION . . 

II.  WAVE  NUMBERS  AND  SCALING . 4 

III.  ACOUSTIC  DOUBLET  FORMATION . S 

IV.  QUADRUPLET  EXPANSION . * 

V.  REFLECTION  COEFFICIENTS  FOR  A  SLOW  BOTTOM . 11 

VI.  RESULTS . 16 

APPENDIX  A . 21 

BIBLIOGRAPHY . 

INITIAL  DISTRIBUTION  LIST . 3o 


IV 


INTRODUCTION 


I  . 

This  thesis  is  a  continuation  in  the  examination  ot  the 
acoustic  pressure  field  in  a  wedge  snaped  ocean.  As  the 
focus  of  naval  operations  shifts  to  the  littoral  regions, 
more  attention  is  being  paid  to  the  problem  ot  A.S.W.  m 
acoustically  shallow  water,  i.e.  regions  where  the  sound 
paths  have  multiple  interactions  with  the  surface  ana  bottom. 
To  date,  the  simplest  propagation  model  has  been  the  metnoa 
of  images . 

In  its  present  form,  this  model  assumes  isospeed  water, 
a  pressure  release  upper  surface,  and  a  lossy,  penetrable 
bottom.  tor  parallel  surtaces,  the  source  and  its  images 
form  a  vertical  array  with  the  actual  source  at  the  center. 
Furthermore,  it  tne  source  is  near  the  suriace  it  can  he 
combined  with  its  negative  (180°  out  ot  phase)  surtace 
reflection  image  to  form  a  dipole  source.  The  column  ot 
images  can  then  be  considered  a  column  ot  doublet  images. 

The  equations  used  in  this  thesis  are  the  far  tiela 
approximations  of  the  acoustic  doublet.  An  unbalanced 
doublet  is  defined  as  having  two  sources  of  nearly  equal 
amplitude  but  opposite  phase  separated  by  a  distance  d.  it 
the  amplitudes  differ  greatly,  or  if  the  phase  difference  is 
not  180°;  then  the  doublet  approximation  is  not  valid. 

When  the  method  of  images  is  used  in  a  wedge  shaped 
ocean,  the  vertical  column  of  doublet  images  becomes  a  circle 
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centered  on  the  apex.  In  previous  work  the  upper  and  lower 
image  doublets  were  summed  individually  with  great  success. 
This  thesis  attempts  to  group  upper  and  lower  douoiets  into 
double  doublets  or  quadruplets. 

This  method  can  only  be  usea  tor  extremely  small  wedge 
angles.  With  larger  angles  the  sound  paths  from  the  upper 
and  lower  doublet  images  will  be  to  dissimilar,  and  the 
quadruplet  approximation  will  be  invalid.  Since  typical 
slopes  for  most  of  the  worlds  continental  shelves  are  about 
3°,  a  good  analogy  would  be  that  ot  using  a  ship  mounted 
active  sonar  searching  for  a  diesel  submarine  in  the  shell 
region . 
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II.  WAVE  NUMBERS  AND  SCALING 
Sound  propagation  in  a  shallow  channel  with  a  pressure 
release  surface  parallel  to  a  rigid  bottom  has  a  propagation 
wave  number  for  each  mode  defined  by  the  following  equation 


ft  1 

H  f  (0) 


Where  H  is  the  depth  of  the  channel . 
is  the  critical  angle  found  from 


U) 

For  a  fast  bottom,  6 


3in0e” 


£i 

^2 


12) 


An  analogous  reference  angle  0£  tor  a  slow  bottom  is 


tan0=— * 
c. 


U) 


If  the  depth  H  is  replaced  by  the  scaling  distance. 


R 


R 

x 


14) 


and  the  bottom  is  tilted  at  angle  fl;  tnen  eq .  1  for  tne 
cutoff  value  ot  the  lowest  mode  with  a  fast  bottom  becomes 


FCx= 


ft 

2sin0tanf) 


and  for  a  slow  bottom 


(5) 


Xxm _ J _ 

2tan0tanfJ 
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III.  ACOUSTIC  DOUBLET  FORMATION 


The  sound  propagation  for  spherical  spreading  of  a  point 
source  is 

(?) 

When  two  point  sources  of  opposite  phase  are  combined  to  form 
a  doublet,  equation  7  can  be  expanded, 

jv  Af  _  A-  ej<*t-*r,) 

r  r. 


Where  the  "+"  and  subscripts  refer  to  the  upper  ano  lower 

sources  respectively.  Tnese  values  can  be  expanded 


i',=r+Ai  A.^A0+J\A 

r.«riix  A_*A0-AA 


i  y  ) 


Where  the  delta  values  indicate  tne  incremental  difference 
tor  each  source  compared  to  a  tneoretical  point  source 
located  exactly  between  them.  The  pressure  equation  can  then 
be  rearranged  to  yield 


1  +  Ad  1-Ad 

[( _ d*-)ejkAr-( - ~)e 

1  i -Al  i+ A£ 

r  r 


uo> 


An  acoustic  source  located  near  the  surface  has  a 
surface  reflection  image  of  equal  amplitude  and  opposite 
phase.  In  the  far  field  the  image  pair  is  treated  as  a 


single  doublet  source.  For  a  balanced  source,  tne  pressure 
equation  can  be  approximated  by 


[e-J*U-Ar)  j 


lll> 


As  seen  in  fig.  2,  Ar  is 


Ar--^sin0 

2 


(12; 


or 

Ar^rjSinYSintf 

(13) 

combining  eqs .  11  and  13,  yields  the  pressure 

P=2j^sin (/nrsinYsir>0)  eJ(wt'*r)  i  ih  > 

To  develop  an  equation  tor  the  unbalanced  doublet,  the 
amplitude  relationships  expressed  in  eq.  9  must  be  rewritten 
as 


0  2  2 


i  io ; 


So  that  for  an  unbalanced  doublet 


2-—  [sin(irr1sin0)  -j^~- cos  (tocsin©)  ]  ci<wt'kr)  i  16; 

1  Ao 


The  distance  from  the  center  ot  the  doublet  to  the 
receiver  can  be  calculated  using  the  law  of  cosines. 


b 


U7  ) 


i*-x*+iI-2z1i2co&  (2n|J±d) 

Where  the  ranges  r2  and  r2  are  measured  from  the  apex  to  tne 
source,  and  from  the  apex  to  the  receiver  respectively.  Tne 
angle  6  is  measured  from  the  surface  down  to  the  receiver. 
The  range  from  each  successive  doublet  m  a  wedge  of  angle  3 
is 


Ub> 


Where  the  angle  m  the  cosine  term  refers  to  the  angle  or 
each  reflected  image;  those  with  a  "+"  are  the  upper  images, 
and  those  with  a  are  the  lower  images.  This  equation  can 

be  rearranged  to  yield 


r»-<r=-r=’N 


Or  r 

~  '  x  ‘  4 


1  + — (1  -COS  (2/jp±6) 
r.-r  2 


( iy ) 


The  first  order  Taylor's  series  expansion  can  tnen  be  usea  to 
simplify  eq .  i 9 


rB*  (r1-r2){l  +  — l-±h-  [i-Cos  (2nPi6)  ]} 

ri-r2 


(20) 
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IV.  QUADRUPLET  EXPANSION 


The  upper  and  lower  image  doublets  can  be  combined  to  torm 
a  quadruplet.  The  complimentary  doublets  are  180°  out  or  phase, 
and  the  pressure  of  the  new  quadruplet  is 

Pn*— ^2(sin[irr1ysin(2r5p+6) J  e'JitAr*  > 

-sin[!c.r1Ysin{2.np-&)  3  e~jlr*r 

In  this  case,  both  of  the  interior  exponents  involving  Ar  are 
negative.  The  Ar  terms  are  derived  trom  a  comparison  to  tne 
range  trom  the  primary  doumet,  thereby  ensuring  that  all  pnase 
angles  are  calculated  with  respect  to  a  common  reterence  point. 
To  find  the  Ar  values  used  m  eg.  21,  an  approximation  or  eg.  iu 
gives 

Ar>r.4 -  (^2_ri)  =— [1-cos  (2nPi5)  ]  12/; 

T  i  ~  r  X 

The  cosine  term  can  be  expanded  and  approximated  ror  small 
values  ot  6. 

cos  (2nP±6)  *cos  (2/iP)  ±6sin(2xiP)  ( 2d; 

Equation  21  can  now  be  expanded  into  the  unbalanced  quadruplet 
equation . 


Prr 


.  2 


- J k— 

gjiwc -kz)  q  ■‘a 


(l-co«  (2ap) 


-jk  -xZ*  dcin  '2 nQ) 

{cin[/tr1Ycin(2np+8)  3  c  r*‘r 1 

Ik  Z^2  5ein'2n(h 

-sin  [/cr:Ysin  (2np-8)  3  <?  Jr*'1  I 


(24; 
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It  is  convenient  to  define 


$=Jcxxn  (i-cos  (2/iP) ) 
dmkT1&8in  (2/ip) 


So  that  eq.  24  becomes 

n  mUel^t-kz)  e-# 

" 

U„lsin  Mcr,ysin(2JiP+6)  ]  -i—coB  l*r,ysin(2np+6)  J  ]  ©'•'■“*  ( 2b> 

/i 

-Bn [sin [icrxY3in  (2nP  -8)  3  -j—£coz  [ic^ysin (2np-8)  3  3  clHd) 

B 

The  values  tor  Aa  and  Ah  must  be  determined  rrom  tne  rer lection 
coefficients . 
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V.  REFLECTION  COEFFICIENTS  FOR  A  SLOW  BOTTOM 


The  reflection  coefficients  used  m  this  thesis  are 
derived  from  the  Rayleigh  reflection  coefficient 

bc-sinOj/sinBj 

bc+sin0t/sin0i 


Where 


b~—  (281 

PA  Cx 

The  angles  are  the  grazing  angle  0t  and  angle  or  transmission 
0t.  Therefore,  6C  is  defined  as 

c-  -arccosO.  (2*1 

Using  trigonometric  rel ationsnips ,  eq.  2/  can  oe 
rewritten , 


bc-J^Tc2 cos2©.  /  6  in0 , 

— i - 1  (JU1 

bc+^/i  -  r^co820i/s  in0a 

In  addition,  the  term  under  the  radical  can  be  rewritten. 


^1  -c^cos^-c^sin^-sin3©,. 


UD 


Combining  eqs .  30  and  31,  and  rearranging  yieias 

W1~tsinee/sinei]a  ^  b-ji-1/x1 
b+ y  1  -  £sitt0c/siiiOi]  2  b+J  1-1/a5 
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Where  x  is 


_  binOj 

X  sin0t 


(3u) 


For  a  slow  bottom  with  a  very  small  angle,  eg.  3/  can  oe 
approximated  by 


Rme 


7^' 

(l-C* 


The  argument  of  the  exponential  is 


r_  2bc  _  2pa/pj 


VT^  yjKCjc^-1 


(34) 


(30) 


So  that  this  approximation  ot  the  Rayleigh  reflection 
coetticient  can  toe  written  m  a  more  concise  torm 

R-e  1164  (3b> 

Since  each  ray  intersects  the  bottom  and  surrace  multiple 
times,  a  product  of  reflection  coefficients  aiong  each 
individual  path  is  required.  This  cumulative  coefficient 
wiii  taxe  the  torm 


2  r.  1 

r„=  n  ^[®i“ (2u-in) pj  (3/) 

•■1,3,... 


and  the  product  yields 


aa®ie-fl3»P 


(  3b  ) 


For  very  small  angles,  this  can  be  approximated  toy 


12 


(39) 


The  angles  of  incidence  are  needed  to  calculate  tne  product. 
The  angles  for  the  bottom  images  are 

0o-2  (ii- 1)  P  (4U) 

Using  the  law  of  sines,  the  angle  at  the  receiver  is  defined 

e=  {iJzq)  sin(2/2p*5)  (41) 

The  angle  of  incidence  of  tne  ray  from  the  center  of  tne 
doublet  image  to  the  receiver  intersecting  the  apparent 
bottom  is 

0=2  (n-1)  (J*6+fc  (4^;; 

However,  the  upper  and  lower  images  m  each  doublet  nave 
slightly  different  angles  of  incidence  with  the  apparent 
bottom 

04=  (2u-l)  (4j) 

The  upper  and  lower  angles  at  the  receiver  tor  each  aouhiet 
are 


*****  Ui/r0)  ycos  (2iip+6)  (44) 

The  individual  angles  of  incidence  are 

(l+  Ui/re)cos  (2nP+fi) )  *  (ij /r0)  ycos  (2ii0)  ( 4b) 

If  the  following  variables  are  used  in  the  coefficients  as 
approximations  of  the  sine  ot  the  angles  of  incidence 
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(  4b; 


an-l*— 1  ccs2n$ 
T0 

ba~  ~  coc2nfi 

ro 

then  eq.  38  becomes 


9Lnr»re^e^  (4  n 

For  each  doublet  image  there  is  an  upper  and  lower  image. 

The  upper  doublet  is  denned  by  AQ  with  tne  upper  image  A. 
and  the  lower  image  At.  Likewise,  the  lower  oouDiet  bD  nas 
the  upper  image  B.  and  tne  lower  image  bt. 


!MB8a 


A.*9t„e  ~a’~e 


B  -*,ena»4aer^a 
£* = Rnt#  "*-**  ts  ~nb-ym 


14b) 


A„  and  Aa  detined  by  eq .  lb  can  be  combined  witn  eq.  4b 


An*K„cosh  (nt>nya )  e  na"#B  fin»Wrcosh  (nbrya )  e 
AA-R^inh  (nbnya)  e'n*°4*  AB*9l„aiuh  {ub„ya )  e 

When  combined,  they  yield 


na-fta 


AA. 


■tanh  ( nbnya ) 


A^*iranh(nbay«) 

Bn 


(i>0) 


The  result  is  the  full  quadruplet  equation 


P-p  +V  acoshb  e*'ut~kr)e'J* 

1  &  z 

[©-•-Ji»dgin( kr1ysin(2n$*6) )  -jtanhb  cos  (ArjYsin  (2nP+4) ) 
-e«*ii»a8in(icxlYsin(2i2p-5) )  -jlanhb  cos  (Ax1ysin(2/ip-&) )  ] 
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LU 

a 


KESULTb 


VI  . 

Three  runs  at  dirterent  ranges  were  maae  to  compare  tne 
output  or  tne  quadruplet  model  with  tnat  or  tne  “UK'itxT1 
program.  For  each  run,  the  source  angle  is  1°,  and  tne  wedge 
angle  is  3°.  The  ratio  ot  tne  sound  speed  m  seawater  to 
that  of  the  bottom  is  1.01695,  and  the  ratio  or  tne  density 
of  seawater  to  that  ot  the  bottom  is  0.7.  This  is  equivalent 
to  a  slow  bottom  consisting  or  a  combination  ot  ciay  and 
ooze;  conditions  very  common  in  continental  shei i  regions. 

The  runs  are  tabulated  at  0.308  increments  ot  receiver  angie 
from  the  surface  downwards.  The  results  are  snown  m  taDies 
1  through  3.  At  a  receiver  angle  of  zero  all  tnree  runs 
result  m  an  amplitude  ot  zero,  as  they  snoulo. 

At  the  shortest  range.  Run  3,  the  results  ot  the 
quadruplet  model  differ  signi £  leant  1 y  trom  those  or  UKTEXT 
The  percentage  oirterence  in  amplitude  is  consistently  aoove 
10%.  With  the  exception  of  tne  3s  receiver  angle,  the 
average  percentage  difference  in  amplitude  tor  tne  otner  two 
runs  is  well  below  10%.  The  tar  neio  doublet  approximation 
improves  at  longer  ranges,  and  the  approximations  tor  Ar 
values  become  more  accurate.  Also,  the  difference  in  phase 
angle  is  much  greater  at  the  shortest  range;  half  the  phase 
angles  in  Run  3  differ  by  greater  than  x/4,  but  in  the  other 
two  runs,  only  the  phase  angie  at  a  3°  receiver  angle  diners 
by  more  than  that  amount. 


Each  ot  these  data  sets  is  accompanied  t>y  sample  piots 
of  the  quadruplet  pressure  amplitude  and  tne  retiection 
coefficient  for  each  quadruplet.  In  each  case,  the  plots 
were  taken  with  a  receiver  angle  of  0.90°.  The  magnitude  ot 
the  complex  quadruplet  pressure  drops  oft  very  quickly  witn 
the  quadruplet  number.  It  is  negligible  by  the  btn 
quadruplet.  This  is  a  result  ot  the  effect  of  tne  rapia 
(exponential)  decay  ot  the  Rayleigh  retiection  coelticient 
with  increased  angle  of  incidence  ot  the  higher  number  ot 
quadrupl et . 

In  practical  terms,  this  means  that  the  quadruplet 
expansion  becomes  less  accurate  with  larger  receiver  and 
wedge  angles.  As  these  angles  increase,  the  ditrerence  in 
the  angles  of  incidence  between  the  upper  doublet  sound  patn 
and  the  lower  doublet  sound  patn  win  increase 
proportionately.  When  these  differences  are  raised  to  the 
exponential  power  in  the  retiection  coetticienr s ,  tne 
differences  in  amplitude  between  upper  ana  lower  douDiets  at 
the  receiver  will  be  quite  large.  However,  the  fundamental 
assumption  in  the  quadruplet  expansion  is  that  the  upper  and 
lower  doublet  images  are  equivalent  in  amplitude,  but 
opposite  in  phase.  Therefore,  the  quadruplet  expansion  is 
only  valid  tor  very  small  angles  and  i ong  ranges. 
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TABLE  1 

BETA5 3 , GAMMA=1 ,R1=1,R2=100, DENSITY  RATIO5 . 9 , SPEED  RAT I 0=1 .  01095 
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“If 

II 

A: 

II 

II 
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AMP. 

1 

1 

AMP. 

PHASE 

|  PHASE. 

II 
•  1 

PHASE 

II 

li 

0.3° 

0.00010 

0.000092 

1 

1 

8 . 4 

1 

|  17 . 6 

1 

| -23 . 9 

II 

41 . 3 
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APPENDIX  A 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%QUADRUPLET  EXPANSION 

%pquad.m 

%MATLAB 

%Quadruplet  expansion  ot  the  acoustic 
%pressure  field  in  a  wedge  shaped  ocean. 

%Michael  Joyce 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

B= input (  wedge  angle  in  degrees*'): 

G=input( 'source  angle  from  surface  in  degrees* ’ ) 
D=input (' receiver  angle  from  surface  in  degrees* 
rho=input( 'water  to  bottom  density  ratio*'); 
cc=input( 'water  to  bottom  sound  speed  ratio=’>; 
rl=input( ’ range  of  source  from  apex*'); 
r2*input (' range  of  receiver  from  apex= ’ ) ; 

Nl=f ix( 90/B) ; 

B=B*pi/ 180 ; 

G=G*pi/180; 

D=D*pi/lB0; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%CALCULATE  SCALING  FACTOR 


2? 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

tb=tan(B) ; 
if  cc<l, 

tl=acos{ cc ) ; 
t2=sin( tl ) ; 

else 

tl=acos(l/cc) ; 
t2=tan( tl ) ; 

end 

kl=2*t2*tb; 
k=pi/ kl ; 

%CALCULATE  CONSTANTS  AND  PRIMARY  DOUBLET 
%%%%%*%%%%**%%%%%%%%*%*%%%%%%%%*% 

al  =  2*(l/rho)/(sqrt(cc"2-l ) )  ; 

r=sqrt ( rl  2+r2  2-2*rl*r2*cos ID)); 

r3=abs(r2-rl) ; 

r4=rl*r2/r3  2; 

mu=r2/abs( r2-rl ) ; 

q=k*rl *G*D; 

pl= ( 2/r )*sin(q) ; 

f =0 ; 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%QUADRUPLET  SUMMATION 
%%%%%%%%%%%%%%%%%%%%%%%%%%*%%%%%% 

for  n=i:l:Nl, 
s(n)=n; 
th(n)=2*n*B; 
d(n)=k*rl*D*sin(th(n)); 
phi (n)  =  K*rl*mu* ( i-cos i th( n ) )  )  ; 

%%%%%%%%%%%%%%%*%%%%% 

^REFLECTION  COEFFICIENTS 

%*%%%%%%%*%%%%%%%%%%%%*%* 

ruln)=sqrt( ri  2  +  r2  2- 2*ri*r 2*cos ( th ( n ) +  1> )) 

r 1 ( n) =sqrt ( rl  2  +  r2' 2-2*rI*r2*cos( t h ( n ) - D ) ) 

rbin)=rl/ruln) ; 

r6(n)=rl/ rl (n;  ; 

a.un)  =  i  +  r:>ui;*cosithin;  ) ; 

bl(n)=r6(n)*cos( thin;); 

rr tn;=-b*ai »(n  2); 

R(n)=-exp( rr (n) ) ; 
a(n)=n*al*D*al(n) ; 
b(n)=n*al*G*bl(n) ; 


%PRESSURE 


ul(n)=cosh(b(n) ) ; 
u2(n)=exp(-D*phi(n) ) ; 
u3(n)=(-2/r)*R(n)*ulln)*u2(n) ; 

%%%%%%%%%%%%%%%%%%%%%%%%% 

%UPPER  IMAGE 

%%%%%%%%%%%%%%%*%%%**%%%% 

vHn)  =  sin(th(n)+D)  ; 
v2(n)=k*tl*G*vHn) ; 
v3(n)-sin( v2(n) ) ; 
v4(n)  =  cos(v2(n) ) ; 
v5(n)=-a(n)-3*mu*d(n) ; 
vb(n>  =  exp(  v5^) ) ; 
v7(n)=tanh(b(n) ) ; 

vb(n)=v6(n)*(v3(n)-j*v>in)*v4(n)) 

%%%%%%%%%%%%%%%%%%%%%%%%* 

% LOWER  IMAGE 

%%%%%%%%%%%%%%%%%%%%%%%%% 

wl(n)=sin(th(n)-D); 
w2(n)=k*rl*G*wl(n) ; 
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w3(n)=sin(w2in) ) ; 


w4 1  n )  =  cos 1 w2 1 n )  ) ; 
w5(n)=a(n)+3*mu*d(nj ; 
w6(n)=exp( w5(n) ) ; 

w7(n)=w6{n)*(w3(n)-;j*v7(n)*w4(n) )  ; 

%%%%%%%%%%%%%%%%%%%%%%%%% 

%SUMMAT ION 

p(n)=u3(  n )  *  (  v8 1  n )  -w7  inn; 
t  =  t  +  p(n )  ; 
pr(n)=real(p(n) ) ; 
pinHn)  =  imag(pm;  ) ; 
pz(n)=sqrt(pr(n)  2+pim(n)  2); 
rrmn  >  =  - 1 x  K  i  n  ;  ; 

end 

piotis.pz;  ;gna; 

xlabel 1 'Quadruplet  Number'); 

ylabel ( 'Magnatude  ot  Quadruplet  Pressure'); 

title( 'Quadruplet  Pressure'); 

pause; 

plot(s, rm) ;grid; 

title! ’ Ray leigh  Reflection  Coet ticient  ) ; 
xiabel ( ’Quadruplet  Number'); 
yiaoeii  Coer t lcient ' ) ; 
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pause ; 


ps  =  pj.  +  t ; 

ph*imag(ps)/reai  (ps ) ; 
phl=atan(ph) ; 
phase- 180 *phl/ pi 
t  =  f  ; 
pl=pl; 
ps=ps : 

P=sqr t ( rea i ( ps )  /+imagips>  /) 
end  ; 


62 


BIBLIOGRAPHY 


Baek  ,  C .  ;The  Acoustic  Pressure  in  a  Hedge  shaped  Hater  Layer 
Overlying  a  Fast  Fluid  Bottom;  Master's  Thesis,  Naval 
Postgraduate  School,  Monterey,  CA,  Marcn  I9b4 

Bradshaw,  N.;  Propagation  or  Sound  in  a  Fast  bottom 
Underlying  a  Wedge  Shaped  Medium;  Master’s  Thesis,  Navax 
Postgraduate  Scnool ,  Monterey,  CA,  Sept.  i960 

Buckingham,  M.J.;  "The  Theory  of  Three-dimensional  Acoustic 
Propagation  m  a  Wedgelike  Ocean  with  a  Penetraoie  Bottom;  ' 

J.  Acoust.  Soc .  Am.,  82(1),  July  1987 

Coppens ,  A.B.,  humpnnes,  M.,  Sanders,  J.V.;  Propagation  oi 
Sound  out  of  a  Fluid  Weage  into  an  Underlying  Fiuia  suostrate 
ot  Greater  Sound  Speed;"  J.  Acoustic.  Soc.  Am., 

76( 5) , November  1984 

Coppens,  A.b.,  rrey,  A.K.,  Mnsier,  l  .  c. .  ,  banaers,  J.v.; 
Fundamentals  or  Acoustics ,  Third  Edition;  John  Wiiey,  ana 
Sons, Inc.;  Monterey,  CA ,  i9«u 

Coppens,  A.B.,  Sanders,  J.V.,  Xoannou,  G.I.,  Kamawamura, 

M . ; Two  Computer  Programs  tor  the  Evaluation  or  the  Acoustic 
Pressure  Amplitude  at  the  Bottom  or  a  Hedge  Shaped  Fluid 
Layer  Overlying  a  Fast,  Fluid  bait -Space;  Technical  Report 
#NPS-61-79-002 , Naval  Postgraduate  School,  Monterey,  CA,  Dec. 
1978 

Jensen,  F.B.,  Kuperman,  W.A.;  'Sound  propagation  in  a  weage 
Shaped  Ocean  with  a  Penetrable  bottom;"  J .  Acoustic.  Soc. 

Am.,  67  is).  May  196U 

Jensen,  F.B.,  Tindie,  C.T.;  'Numerical  Modeling  Kesuits  tot 
Mode  Propagation  m  a  Wedge;"  J .  Acoustic.  Soc.  Am.,  b2vi), 
July  1987 

Kaswandi ,  C.;  A  Computerized  investigation  Using  the  Metnoa 
of  Images  to  Predict  the  Sound  Field  in  a  Fluid  Hedge 
Overlying  a  Slow  Halt-space;  Master's  Thesis,  Navai 
Postgraduate  School,  Monterey,  CA,  Dec.  198/ 

Kim  Jong  Rok;  Comparison  tor  Sound  Pressure  in  a  Weage  Shaped 
Ocean  as  Predicted  by  an  Image  Method  and  a  PE  Model; 

Master's  Thesis,  Naval  Postgraauate  School,  Monterey,  , 

Dec. 1990 


Livmgood,  D.M.  ;  Extension  or  tne  Analytical  Approximation  to 
the  Transmission  of  Sound  m  Shallow  Mater  Using  tne  Image 
Model;  Master's  Tftesis,  Naval  Postgraduate  scnool ,  Monterey, 
CA,  Sept.  1992 

Math  Worxs ,  inc. ;  me  Student  taition  ot  MArLAb;  crentice 
Hail,  Inc.,  Englewood  dirts,  New  Jersey,  u?j62 

Nassopoulis,  G.  ;  Study  or  Sound  Propagation  in  a  heage  snapea 
Ocean  and  Comparison  witn  other  Methods;  Master's  rnesis, 
Naval  Poastgraduate  Softool  ,  Monterey,  ca  ,  June  iw2 

Paliatsos,  D.;  Computer  Studies  ot  Sound  Propagation  m  a 
Wedge  Shaped  Ocean  with  Penetrable  bottom ;  Master's 
Thesis,  Naval  Postgraduate  School ,  Monterey,  CA,  Marcn  19«y 

Personal  Communications  with  A.b.  coppens  and  J.V.  Sanders 
Naval  Postgraduate  School,  Monterey,  CA,  June-Sept.  19*.j 


DISTRIBUTION  Li S*1 


1.  Detense  Technical  Intormation  Center 
Cameron  Station 

Alexandria,  VA,  22304-6145 

2.  Library,  Code  0142 

Naval  Postgraduate  School 
Monterey,  CA ,  93943-6002 

3.  Dr.  A . B .  Coppens,  Coae  PH/CZ 
Department  ot  Physics 

Naval  Postgraduate  School 
Monterey,  CA,  *6*4j-suQ2 

4.  Lt  Michael  D.  Joyce 
ASWOC  Sigoneiia 
PSC  812,  BOX  32 /0 
FPO  Ah  0yo^/-62'7u 


35 


